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ABSTRACT 

We investigate the reheating of the very late thermal pulse (VLTP) object V4334 Sgr (Sakurai's Object) using radio observations 
from the Very Large Array, and optical spectra obtained with the Very Large Telescope. We find a sudden rise of the radio flux at 

5 and 8 GHz — from < 90 /jJy and 80 ± 30 /jJy in February 2005 to 320 /jJy and 280 /jJy in lune 2006. Optical line emission is 
also evolving, but the emission lines are fading. The optical line emission and early radio flux are attributed to a fast shock (and not 

' photoionization as was reported earlier) which occurred around 1998. The fading is due to post-shock cooling and recombination. 

The recent rapid increase in radio flux is evidence for the onset of photoionization of carbon starting around 2005. The current results 
1 1 . indicate an increase in the stellar temperature to 12 kK in 2006. The mass ejected in the VLTP eruption is M c j > 10~ 4 M Q , but could 

be as high as 10~ 2 M Q , depending mainly on the distance and the clumping factor of the outflow. We derive a distance between 1.8 
and 5 kpc. A high mass loss could expose the helium layer and yield abundances compatible with those of [WC] and PG1 159 stars. 

f-^. ' Key words. Stars: individual: V4334 Sgr - Stars: mass loss - Stars: evolution - Planetary nebulae: general - dust, extinction 
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the VLTP, during which the star became hydrogen poor and 
carbon rich. The speed of evolution from the VLTP to max- 
imum brightness was surprising as pre-Sakurai models pre- 
dicted timescales of > 100 yr. The much faster evolution has 
been explained as c aused by suppressed convec ti on under ex- 
plosiv e conditio ns dLawlor & MacDonaldl 120031; lHaiduk et alj 
120051) . However. [Miller Bertolami et al.l d2006l) find that the nu- 
merical time step is important, and that smaller time steps are 
sufficient to predict much fa ster evolution . The p roblem of con- 
vection is studied further bv lHerwig et al.l ((2006). The reheating 
timescale after maximum brightness is an important constraint 
to test the models. This evolution can be tra ced best by observ - 
ing the ion i zation of the surrounding n ebula. lEvres et al.l dl999). 
iTvne et all d2000l) . and iKerber et all (|2002|) detected emission 
from atomic and ionized species around V4334 Sgr which is 
evidence for ionization . The existence of ionization was con- 
firmed by lHaiduk et al.1 d2005l) based on radio observations. In 
this Letter, we present new radio and optical observations, show- 
ing that this early ionization was due to a fast shock. Evidence 
for the onset of photoionization of carbon is also found, indi- 
cating that this started around 2005. We present photoionization 
models of the ejecta and compare the results to various VLTP 
evolutionary tracks. 



2. Observations 



1. Introduction 

\Q • 

■ Helium shell flashes dominate many aspects of the evolution of 
' AGB (Asymptotic Giant Branch) and post-AGB stars. Helium 

f^) , flashes (or thermal pulses) are common on the AGB but the ef- 
1| ■ fects on the surface properties are mitigated by the thick hydro- 
. £h ! gen envelope. Very late thermal pulses (VLTPs), which can oc- 

■ cur after hydrogen burning has ceased, lead to large and rapid 
$— i ' changes in the star and are the best way to constrain the physics 
5^ , of these thermonuclear eruptions. The VLTP in Sakurai's Object 

(V4334 Sgr) has provided the first opportunity to observe such 
an event with modern instrumentation, and allows us to study 
poorly understood aspect s, such as the post- V LTP mass loss and 
the mixing length theory (Haiduk et al. 2005[). 

Very few VLTP event s have been ob served: only V4334 
Sgr (iPuerbec ket al.lll997|) and V605 Aql dClavton et al.ll2006t 
iLechner & Kimeswengerl |2004|) were discovered during their 
high-luminosity phase (respectively in 1996 and in 1918). CK 
Vul (i n 1670) is suspected to represent a third case (Haidu k et al.l 
|2007|) . FG Sge also shows eviden ce for current post- VLTP evo - 
lution but its status is uncertain dJefferv & Schonbernetj|2006l) . 
Another five central stars of planetary nebulae (PNe) show ev- 
idence for historical VLTP eruptions, b ased on the pr esence of 
hydrogen-poor gas near the central star (Zijlstra 2002). 

V4334 Sgr showed a high luminosity phase with a cool 
stellar atmosphere (r e jf ~ 6000 K) within a few years after 



* Based on observations collected at the European Southern 
Observatory, Chile (programmes 71.D-0396, 75.D-0471, 77.D-0394). 



Radio emission We have been monitoring the radio flux 
of V4334 Sgr to test the temperature evolution of the star. 
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Observations were carried out with the VLA array between 2004 
and 2006, at 5 GHz and 8 GHz. We observed in the hybrid BnC 
(2004) and AnB (2005, 2006) arrays. This reduced confusion 
with emission from the old, extended PN. Integration times were 
3-5 hours, with 1331+305 as the flux calibrator and 1733- 130 as 
the phase calibrator. The assumed flux density of 1331+305 was 
4.66 Jy at 5 GHz. Flux calibration was interpolated to V4334 
Sgr. Observational parameters are summarized in Table 1 . The 
resultant visibilities were Fourier transformed to the image plane 
and the dirty beam deconvolved. Convolution with a gaussian 
beam fitted to the dirty beam size gave images from which fluxes 
were measured. Natural weighting was used for the imaging. 
The FWHM of the gaussian beams are listed in Table 1 . 

Optical spectroscopy Line emission from the central ejecta in 
the optical wavelen gth regime was first detected in 2001 by 
iKerber et all d2002l) using the FORS1 instrument on the VLT. 
We continued to monitor the object with FORS 1/2 in subsequent 
years. Details are given in Table [2] We averaged observations 
taken within the same semester. The spectra were manually re- 
duced in the standard way using iraf. Sky lines were removed 
using the part of the long slit spectrum outside the old PN. The 
old PN also contributes line emission: this is strongly depen- 
dent on position on the slit, and the old PN lines were removed 
by interpolation of the area around the central source. The re- 
sponse correction was done using various white dwarf flux stan- 
dards. The line fluxes derived from the VLT spectra are listed in 
Table [2] The data we re not corrected for reddening. The 2001 
values are taken from lKerber et alj d2002l) . 



Table 1. Observational parameters and results for the VLA ob- 
servations. PA gives the position angle of the beam. 



Date 


Freq. 


FWHM Beam, PA 


cr 


flux 




[GHz] 




[/<Jy/b] 




5 Feb 2004 


8 


1. 30" x 1.14", +48° 


10 


100 ± 30 


4 Feb 2005 


5 


1.37"x0.96", +68° 


16 


< 90 


6 Feb 2005 


8 


0.78"x0.54", +76° 


14 


70 ±40 


11 Jun2006 


8 


0.73"x0.65", +46° 


17 


280 ± 50 


12 Jun2006 


5 


1. 43" x 1.05", +2° 


19 


320 ± 60 



Table 2. Line fluxes in units of 1 0" 17 e rg cm' 2 s" 1 . The 2001 
fluxes are taken from lKerber et al.l d2002l) . OR means outside the 
observable range. Where appropriate 3-cr upper limits are given. 



Line 


2001.44" 


2003.50" 


2005. 48 r 


2006.29" 


[N i] 5 198, 5200 


8. 


3.2 + 0.5 


2.6 ± 0.5 


2.2 ± 0.7 


[Nil] 5755 


2. 


2.1 ±0.7 


< 1.5 


<2.2 


[Oi] 6300 


30. 


10. ± 1.5 


5.5 ±0.8 


3.9 ±0.8 


[Oi] 6364 


8. 


4. ± 1.2 


2.5 ± 0.6 


2.1 ±0.7 


[Nil] 6548 


40. 


15. ± 1.9 


7.6 ±0.9 


3.7 ± 0.6 


[Nil] 6583 


121. 


42. ±4.2 


21. ±2.2 


13.8 ± 1.5 


[Sn] 6716 


2.5 


< 1.9 


< 1.6 


< 1.3 


[Sn] 6731 


2. 


< 1.9 


< 1.6 


< 1.3 


[On] 7319, 7330 


62. 


20. ±2.1 


6.7 ± 1.6 


4.2 ± 1.0 


[Ci] 8727 


OR 


13. ± 1.5 


<3.9 


4.7 ± 0.8 


[Ci] 9824, 9850 


OR 


150. ± 15. 


OR 


OR 



" 300V+GG375, 4500 - 8000 A. >' 1501, 4900 - 10075 A; 150I+OG590, 6200 
- 10075 A; 1200R+GG435, 5875 - 7365 A. c 300V+GG435, 4800 - 8900 A. 
d 600V+GG435, 4520 - 6850 A; 600I+OG590, 6500 - 9200 A. 



3. Discussion 

3.1. Extinction distance 

The FORS long slit spectra were used to extract the spectrum 
of the old PN, taking care not to include any field stars or the 
recent ejecta. The Ha/H/3 ratio gives the interstellar extinction. 
Assuming a case B ratio of 2. 85, we derive E (B - V) = 0.86, 
somewhat higher than found b v lPollaccol d 1 9991) . The extinction- 
distance diagram for stars within a few arcminu tes from V4334 
Sgr is derived by Kimeswenger & Kerberldl998l) : it shows a lin- 
ear rise with distance up to 1.8 kpc, followed by a constant ex- 
tinction of E(B - V) = 0.9 ± 0.09 between 1.8 and 5 kpc. 

This suggests that V4334 Sgr is lo cated beyond 1 .8 k pc. The 
previously used value of 1.9 kpc (e.g.lEvans et al.ll2006l) can be 
considered a lower limit. iDuerbeck et aL ( 20001) suggest d ~ 3- 
5.4 kpc based on the stellar luminosity. The line of sight reaches 
the scale height of the old disk at d ~ 4 kpc and this would 
provide a plausible distance. 

3.2. Shock ionization 

The line fluxes in Table|2]show a strong exponential decline with 
time (see also Fig. [TJ. The level of excitation is also decreas- 
ing, as the high excitation [O n] lines are decreasing faster than 
lines of lower excitation. The e-folding times are 1.81+0.13 yr 
for [On] (the highest excitation line), 2.29+0.14 yr for [Nn], 
2.4+0.2 yr for [Oi], and 3.7+0.9 yr for [Ni]. This is clearly in- 
consi stent with photoionization by an incr easingly hot central 
star dKerber et al.ll2002l:lHaiduk et al.ll2005l) . The decrease in in- 
tensity as well as the level of excitation is consistent with a sin- 
gle shock that occurred some where before 2 001 and then started 
cooling and recombining (cf. Kafatos 1973). The [N n] line ratio 



in 2001 indicates an electron temperature T e = 3200 - 5500 K, 
depending on the electron density (Table 0. The [Ci] ratio in 
2003 indicates T e = 2300 - 4300 K. These low values indicate 
that cooling had already commenced by 2001. 

The origin of the shock is not clear. No emission lines 
were seen during 1997. The earliest evidence for shock ioniza- 
tion is based on the He i recombination l ine, and dates from 
1998 dEvres et al.l [l999t iTvne etakl l2000b . High-velocity CO 
was seen starting late 1999 or early 2000 dEvres et al .1 12004). 
with an outflow velocity around 290 km s _1 . Optical shock emis- 
sion with velocities up to -350 / +200 kms~' was first seen in 
2001 dKerber et al.| [2002). The s hock occurred close t o the time 
of the onset of dust formation dDuerbeck et al.ll2000l) . Neutral 
molecule chemistry is slow in a hydrogen-defic ient environ- 
ment, while ion-molecule chemistry is much faster (Wood s et al.l 
120021) . Hence UV radiation from the shock may have triggered 
dust formation elsewhere in the nebula through ionization of key 
molecular species. If the ejecta are bipolar, the dust formation 
could be taking place in the dense equatorial region (possibly 
an equato rial disk), while the shocks are internal to a fast polar 
outflow dEvans et al.ll2006l) . 

3.3. Photoionization 

The radio flux was constant or decreased slowly between 2004 
and 2005. This is consistent with the model of an early fast 
shock: the radio data trace the post-shock cooling and recom- 
bination of the shock and not a photoionized region. However, 
the 2006 observations show a strong increase in radio flux, by 
a factor of four within 16 months. The only optical line which 
appears to have strengthened is the [C i] 8727 line, but this is 
uncertain due to the low S/N of the 2005 data. The increase in 
radio flux in 2006 is therefore attributed to newly photoionized 
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2002 2003 2004 2005 2006 
epoch 

Fig. 1. Evolution of emission line strength. The lines show the 
least-squares fit to the data. 



carbon which has the lowest ionization potential of all the ma- 
jor constituents in the gas. This process would not produce any 
noticeable changes in the optical or mid-IR spectrum. 

We ran Cloudy models to investigate the origin of the ra- 
dio flux using versi on c07.04.01 of the code, last described by 
Ferla nd et alJ (Il998h. The basi c models are similar to those de- 
scribed in lHaiduk et alj d2005l) . with the following c hanges. (1) 
The an gula r diameter of the du st shell is discussed in lTvne et alj 
d2002l) and lEvans etaf] d2006l) . Extrapolating their data to 2006 
leads to an angular diameter of ~ 100 mas. However, with such 
a small diameter the plasma would become optically thick at 
5 GHz before the necessary emission measure could be reached. 
Based on the FORS2 1200R spectral image taken in 2003 (see- 
ing 0.54") we measured a deconvolved (gaussian) FWHM of 
0.27" for the [N n] emitting regio n. This indicate s that the true 
diameter must be 0.3 - 0.5 arcsec (van Hoof 2000). We modeled 



both l imits of this range. (2) Abundances from lAsplund et alj 
( 1999, October 1996 values) were used, ex cept for carbon where 
we used C/He = 0.1 by number (the value lAsplund et all (l999) 
used for their model at mospheres). (3) For t he ste llar spectrum 
we used the models of ICastelli & Kuruczl d2004l) . (4) We as- 
sumed a clumpy medium with a filling factor / of 0.1 or 0.01, 
and a 1/r 2 density distribution. (5) The total dust mass was the 
same in all models with the same distance (see Sect. 13.5b . We 
list four models in Table [3] The implied gas-to-dust mass ra- 
tio is ~ 580 V7 for the d = 1.9 kpc models and ~ 870 V7 for 
the d — 4 kpc models. The ionization of carbon begins around 
r e ff =10 kK. The radio flux (without the shock emission) as a 
function of T e ff is listed in Table [4] The observed radio flux is 
consistent with an increase in T e ff from ^11 kK in 2005 to ~ 
12kKin 2006. 

3.4. Stellar evolution models 

Three evolutionary tracks have been published which can 
fit the rapid post He-flash evolution. iHerwigl d200ll) and 
Haiduketal. (2005 ) proposed sup pressed convection after the 
flash. lLawlor & MacDonaldl (120031) showed that with this as- 
sumption, the evolution could be closely matched; they first pre- 
dictedthe 'double loop' evolution. Mill er Bertolami et alj f2006) 
find that the initial fast changes can also be found in models with 
standard convection, by choosing a very small time step in the 
simulations. This raises the issue of the accuracy of the previ- 



Table 3. Cloudy model parameters. The main constituents of the 
gas are helium (He/H = 250 by number) and carbon (C/H = 25 
by number). The electron density in the ionized region is approx- 
imately 25 «h- Ma and M e j are the dust mass and the total ejecta 
mass. 





model 1 


model 2 


model 3 


model 4 


d [kpc] 


1.9 


1.9 


4.0 


4.0 


L[U,] 


2770 


2770 


12280 


12280 


& [arcsec] 


0.3 


0.5 


0.3 


0.5 


log(r in [cm]) 


14.931 


15.153 


15.254 


15.476 


log(r out [cm]) 


15.630 


15.852 


15.953 


16.175 


log(n H [cm- 3 ]/V/) 


4.45 


3.85 


4.3 


3.7 


M d [10- 6 Mo] 


1.95 


1.95 


8.65 


8.65 


M ej /V/[10- 3 M G ] 


1.05 


1.22 


6.92 


8.05 


Table 4. Predicted radio flux (not including the shock emissio 


at 8.435 GHz in yuJy (valid both for / = 


0.1 and 0.01). 


T e ff model 1 


model 2 


model 3 


model 4 




10000 4 


4 


4 


3 




11000 24 


24 


24 


24 




12000 340 


330 


350 


340 




13000 830 


790 


850 


820 





ous models. In all cases the initial fast evolution is caused by a 
hydrogen ingestion flash at the top of the he lium lay er. 

Fig. [2] shows the double-loop model of Herwigj d200ll) and 
Haiduketal. (2005), with the top inset showing the tempera- 
ture evolution during the first return. The middle inset shows 
the model of Miller Bertolami et al. (2006) and the lower inset 
shows a representative model of Lawlor & Mac Donaldl {2003) 
(not necessarily t heir best-fitting one). The observational dat a 
point s come fromlDuerbeck et all d 1997b . lAsplund etaTI dl999l) . 
and iPavlenko & Geb alle (2003). The last two points are from 
this paper. For broadband photometric data (including pre- 
discovery photometry), we have estimated temperatures from 
bolometric corrections. These corrections are based on normal 
giant photospheres and are therefore uncertain. The low temper- 
ature betwee n 1995 and 2000 is cau sed by an expanding pseudo- 
photosphere dDuerbeck et al.ll2000l) . The underlying star likely 
remained hotter. Such an expanding photosphere is difficult to 
model with a stellar evolution code and predictions for T e g in 
this phase are therefore uncertain. However, this does affect the 
overall reliability of the track. 

The three models differ in the predictions for the reheat- 
ing timescale, but neither fully fits the current observations. 
Herwig's model reheats faster, and the other two models con- 
siderably slower than the observations indicate. The reheating 
timescale in Herwig's model can be adjusted by tuning the core 
mass or the m ixing efficiency in the He-flash convection zone 
dHerwigl20 01). There is a degeneracy between these two param- 
eters and the effects of a smaller core mass and reduced mix- 
ing efficiency are to some degree interchangeable. The mass loss 



during the cool phas e may also be an import ant missing ingredi- 
ent in these models (Iben & Mac Donaldl 19951) . 



3.5. Ejecta mass 

The results from the Cloudy modeling can be summarized in the 
following formula for the ejected mass: 



M ej = 6.0 10~ 



J_ 

0.01 



0.5 



4 kpc 



2.5 





OS 



0.3 



m(He) 
2.5 



Ma, 



(1) 
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where © is the diameter of the nebula and m(He) is the He/C 
mass ratio. The dominant uncertainty in the total ejecta mass 
comes from the distance, the volume filling factor, and the He/C 
mass ratio. In the Herwig models, the He/C ratio quickly evolves 
to a value of 2.5 after the star has reached the red end of the first 
loop. We take this value to determine the minimum ejecta mass, 
required by the photoionization models. Hence a very conserva- 
tive lower limit to the ejecta mass is 6 10~ 4 M at d — 4 kpc and 
10~ 4 M© at d = 1.9 kpc assuming extreme clumping (/ = 0.01). 
For a uniform outflow the lower limits would be a factor of 
10 higher. The corresponding mass-loss rate over a period of 
roughly a decade is very high. It could be comparable to, or even 
exceed, the mass-loss rate during the superwind phase. A mini- 
mum ejecta mass of 5 10~ 3 M is sufficient to expos e the inter- 
shell region (e.g. Fig. 3 in llben & MacDonaldl fl995) where the 
C/He ratio reaches its peak value of 2.66 by mass. This would 
leave a central star with abundances compatible with those of 
[WC] and PG1 159 stars. The linear momentum in the ejecta far 
exceeds that in the stellar radiation (even when considering mul- 
tiple scattering), so that the ejecta must be energy driven. This 
suggests that the ejecta may have been lost instantaneously. 

The dust mass determined from the continuum emission 
has been steadily increasing, starting in 1998 or early 1999 
dTvne et alj|2002l) . In June 2003, the total dust mass indicated 
by the sub-mm continuum was 1.9 x 10~ 6 M at d = 4 kpc. 
Between January and June 2003 the dust mass increased by a 
factor of 1.7. The increase was used to deri ve a mass-loss rat e 
of 1 .4 x 10~ 4 M yr 1 at d = 4 kpc by 2003 dEvans et alj|2004b . 
However, the concept of converting the dust growth rate to a 
mass-loss rate assuming a constant and high dust-to-gas ratio 
(1/75 by mass) can be questioned. The low observed dust tem- 
perature, Td x 400 K in 2003 seems difficult to reconcile with 
the premise that half the dust was newly formed and therefore 
located close to the dust condensation radius. Instead, the in- 
crease in dust mass can represent a continuous dust fo rmation 
or growth in a circumstellar disk (cf. Haiduk et al. 2007). Based 
on the flux in the Spitzer spectrum dEvans et al. 1 120061 we find 
a current dust mass Mj * 8.7 x 10~ 6 M Q at d = 4 kpc, us- 
ing opacities for graphite. The dust mass has increased five- 
fold between June 2003 and April 2005 (when the Spitzer data 
were taken). The average dust growth rate during this time was 
M rf * 3.7 x 10~ 6 M yr" 1 , a factor of 2 higher than in 2003. 

4. Conclusions 

The photoionization of carbon in the ejecta of V4334 Sgr has 
started around 2005. This is evident from a 4-fold increase in 
the 8 GHz radio flux observed in June 2006. Simultaneously 
the optical spectrum shows a continuing exponential decline, 
both in excitation and ionization. This indicates that the emis- 
sion lines are formed in a shock which occurred around 1998 
and is currently cooling and recombining. Based on the Balmer 
decrement in the spectrum of the old PN, we could determine 
that the distance to V4334 Sgr is in the range 1.8-5 kpc, with a 
preferred value ~ 4 kpc. New Cloudy models of the photoionized 
region indicate that the total shell mass is at least 10~ 4 M Q , but 
could be as high as 10~ 2 M (assuming / = 1 and d — 5 kpc), 
depending mainly on the assumed distance and clumping fac- 
tor in the outflow. The models indicate that the star was at a 
temperature of 12 kK in 2006 and is currently evolving with 
a speed ^ 1 kKyr -1 . This confirms that V4334 Sgr not only 
displays a fast return to minimum temperature, but also a fast 
reheating. A comparison of the temperature evolution of the 
central star with observations shows that Herwig's model re- 



i — 1 — 

1999.9 




5 4.5 4 3.5 



log T [K] 

Fig. 2. Comparison of predicted and observed evolutionary 
timescales. 

heats faster dHerwigl l200lt lHaiduketal1l2005l) . and the other 
two models considerably slower than the observations ind icate 
dMiller Bertolami et al . 2006i lLawlor & MacDon ald 2003). The 
dust-to-gas mass ratio may be lower than hitherto assumed, de- 
pending mainly on the clumping factor of the gas. Our analysis 
indicates that dust formation and/or growth is slow and still on- 
going. 
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